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Phase synchronization between two superradiant lasers 


Joshua M. Weiner, Kevin C. Cox, Justin G. Bohnet, and James K. Thompson 
JILA, NIST and Department of Physics, University of Colorado, Boulder, Colorado 80309-0440, 

(Dated: March 24, 2015) 

We experimentally demonstrate synchronization between two distinct ensembles of cold atoms undergoing 
steady state superradiance within a single longitudinal and transverse mode of the same optical cavity. The 
synchronization process is studied first in terms of the time dynamics of re-synchronization when the phase 
alignment of the two oscillators is abruptly broken. We also observe the steady state behavior of the lasers 
as their relative frequency is continuously varied. This system has the potential to realize a non-equilibrium 
quantum phase transition and could inform future implementations of milliHertz linewidth lasers. 

PACS numbers: 42.55.Ye, 32.80.Qk, 37.30.+i, 42.50.Ct 


Phase synchronization of oscillators is a ubiquitous phe¬ 
nomenon, occurring in physical, chemical, biological, and 
social systems fT]. Recent demonstrations of synchronized 
oscillators at the nano-scale have been reported in mechani¬ 
cal 0, opto-mechanical EH, spintronic EE, and electro¬ 
mechanical 13 systems. Synchronization dynamics have also 
been observed in state-of-the-art frequency combs m Ill- 
Synchronization has been extensively explored with classical 
models iioiini, and more recently systems have been identi¬ 
fied in which quantum noise could contribute to synchroniza¬ 
tion dynamics lfT2m4ll . 

Synchronized open quantum systems are of fundamental 
interest in exploring synchronization models 1141 , associative 
memories ca, and quantum computing El. Quantum ef¬ 
fects on synchronization are expected in opto-mechanical ifTSl 
ITtUT^ , optical i20l , and cold atom lfT^IT4ll^ systems. Here 
we study synchronization dynamics in the system introduced 
by Ref. Ql: two laser-cooled atomic ensembles that interact 
through a common cavity mode via steady state superradiant 
emission. 

The synchronization mechanism of steady state superradi¬ 
ant lasers may enable milliHertz linewidth optical lasers that 
are highly insensitive to both technical and thermal mirror 
vibrations Eaiiii, an optical analog of the microwave hy¬ 
drogen maser 12^ . In such a laser, cavity-mediated interac¬ 
tions combined with repumping-induced dissipation cause the 
spontaneous synchronization of the phases of the radiating op¬ 
tical dipoles of individual atoms. In the absence of synchro¬ 
nization, the optical dipoles would quickly dephase due to 
both homogeneous and inhomogeneous broadening, leading 
to weaker incoherent light emission with a linewidth directly 
reflecting the width of the broadened atomic transition. Such 
a narrow optical frequency reference would find a broad range 
of applications in timekeeping, long-baseline optical interfer¬ 
ometry, and precision measurement ll25l . 

In this Letter, we study the synchronization of two distinct 
sub-ensembles of atoms whose relative optical dipole phases 
can be externally controlled. This allows us to abmptly break 
the phase alignment between the two collective optical dipoles 
and watch as synchronization heals the relative phase error. In 
a second set of experiments, we introduce a continuous source 
of phase errors (i.e., a frequency offset) between the two en¬ 


sembles and observe how the competition between the heal¬ 
ing rate and the phase error rate leads to a threshold for syn¬ 
chronization to occur. The experimental results are in good 
agreement with a mean-field model that does not treat quan¬ 
tum noise. 

To form the superradiant laser gain medium, we prepare 
N = 1.2 X 10® ®^Rb atoms at 20 pK within a ID op¬ 
tical lattice in a high-finesse optical cavity with power de¬ 
cay rate k = 27r x 12 MHz and single-atom coopera- 
tivity C = 5 X 10“®. The lasing transition is a Ra¬ 
man transition from the |t) = \5‘^Si/2,F = 2,mF = O) to 
14,) = \5‘^Si/ 2,F = l,mF = O) ground hyperfine states. In a 
dressed-state picture, the effective atomic transition frequency 
is the frequency of the spontaneously emitted Raman pho¬ 
ton. The transition frequency is controlled by Raman dressing 
lasers applied transverse to the cavity axis and tuned 1.3 GHz 
blue of the |t) to |i) = , F = 2) transition. Repump¬ 

ing from 14,) back to |t) at single-atom rate W is achieved by 
applying additional lasers transverse to the cavity mode. The 
repump lasers are not phase matched with the Raman dress¬ 
ing lasers and their spontaneously scattered photons are not 
resonant with a cavity mode. 



FIG. 1. (Color online) Experimental diagram and Raman lasing en¬ 
ergy levels, (a) Two spatially distinct beams (red, blue) dress an 
ensemble of laser-cooled atoms inside an optical cavity, defining the 
two superradiant ensembles a and b. Repumping beams (green) are 
also applied transverse to the cavity, (b) Dressing beams induce Ra¬ 
man decay from |t) to |^). Both emitted photon frequencies (wavy 
lines) are within the linewidth «: of a single cavity mode. The re¬ 
pumping laser returns atoms back to |t) via single-particle repump¬ 
ing at rate W. 
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To create two spatially separate ensembles with indepen¬ 
dently controlled optical dipoles, we apply two Raman dress¬ 
ing lasers that address either the upper or lower portions of the 
total trapped atomic ensemble (Fig. [T]|. This provides inde¬ 
pendent control of the dressing laser phases aa,b, angular fre¬ 
quencies U!a,b, and intensities as parameterized by a resonant- 
Rabi flopping angular frequency fla^b for the |t) to |i) transi¬ 
tion. We can independently set the single-atom Raman decay 
rates ^a,b{~ Stt x 250 Hz) by controlling each laser’s inten¬ 
sity. The relative number of atoms Na^b in each ensemble can 
be controlled by translating the spatial boundary between the 
dressing lasers along the cavity axis. 

Because we utilize Raman transitions for the lasing pro¬ 
cess, the relevant total optical dipole phases that synchronize 
are given by (j)a^b = Va,b + ctafi- Here r]a,b is the phase as¬ 
sociated with the coherence that develops between the ground 
states It) and |t) in each ensemble. Since the dressing phases 
are externally controlled parameters, the cavity-mediated in¬ 
teractions drive changes in the ground state coherences rja^b to 
synchronize the optical dipole phases (j)a,b- 

We first study the dynamics of phase synchronization in the 
time domain for two ensembles with degenerate frequencies 
S = (jJb — uja = 0- The dressing and repumping lasers are all 
turned on for 0.1 ms, during which time the two ensembles 
reach a steady state in which they emit at the same frequency 
and act as a single synchronized superradiant ensemble with 
'Pa = </'&• An electro-optic crystal is used to quickly jump the 
phase ab of the b dressing laser by an amount Aab in 30 ns. 
The timescale of the jump is much faster than the time dynam¬ 
ics of the resynchronization process and effectively creates 
an instantaneous error in the alignment of the optical phases 
4^b ~ 4^a “b Aoib- 

To observe how this phase error heals in time, we allow the 
system to dynamically evolve for a variable amount of time 
Tevoi = 0 to 1.5 /rs before we rapidly extinguish the other 
dressing laser, —> 0. Subsequently, only ensemble b radi¬ 

ates into the cavity mode. We infer the change in 0;, from the 
difference in the phases At/j of the emitted light just before 
the phase jump and just after Tevoi- 

There are several other physical mechanisms that also af¬ 
fect the observed value Aij} that are not directly related to 
the synchronization between the optical dipoles. The pri¬ 
mary contribution to this background phase shift is popula¬ 
tion inversion-dependent cavity frequency pulling 1^ . To 
remove these less interesting contributions, we measure the 
light phase difference A'tjj± for equal magnitude but opposite 
sign phase jumps ±Aat,. The computed differential quantity 
Aijj = — Aijj-)/2 is insensitive to these background 

systematic errors. 

The measured quantity A'0 as a function of the evolution 
time Tevoi is shown in Fig. |^b). Here the phase jump is 
Aab = 90°, and we see that A'0 is also 90° near Tevoi = 0. 
The phase A'l/i then relaxes back toward 0°, settling at an inter¬ 
mediate value such that (j)a = (j)b- The timescale for relaxation 
is close to the repumping rate W~^, i.e., the characteristic rate 


at which phase errors are erased, as discussed below. 



FIG. 2. (Color) Healing of an instantaneous phase error between op¬ 
tical dipoles, (a) Timing diagram and visualization of atomic Bloch 
vectors. Before time t — 0 the two dipoles interact and synchronize. 
At t: = 0, dressing laser phase ai, is jumped by 90°. The ensem¬ 
bles’ interaction begins to heal the relative phase error. At t = Tevoi, 
dressing laser a is turned off (fl^ —>■ 0) so that only ensemble b radi¬ 
ates into the cavity. The difference A-i/i in the phases of the radiated 
light in the grey windows before t = 0 and after t = Tevoi indi¬ 
cates the change in the optical dipole phase Acpi, = Atp. The upper 
panels provide cartoon visualizations of phasors representing the ra¬ 
diated fields (red for a, blue for b, purple for the sum) and Bloch 
vectors, (b) Light phase change A-i/i vs. evolution time Tevoi- The 
solid and open points correspond to experiments with dipole ratios 
Rd = (1.5,4.0) respectively. Vertical solid and dashed lines show 
the characteristic time scale of the respective single-atom repumping 
rates for the two data sets W~^ = (0.77,1.6) fds corresponding to 
(solid, open) data. The solid and dashed curves are simulations for 
the respective data (red for ensemble a, blue for b). 

The equilibrium phase at large Tevoi is mostly determined 
by the ratio of the relative magnitudes of the optical dipoles 
of the two ensembles just before the evolution period. The 
magnitude of each collective dipole is proportional to the 
number of participating synchronized atoms (Na,b) and the 
emitted electric field per atom (c>c y/ja,b)- The relative 
dipole magnitude is then roughly characterized by = 
{NbyAfb)/ {Na's/^) = 1.5 and 4.0 for the solid and open 
data sets in Fig. [^b). A simple model for Tevoi W~^ and 
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Aa = 90° predicts that A'ljj will relax to A'ljje = ta,n~^{Rd)- 
The steady state phase given by the simulation is A-ipn- For 
the data with more balanced populations (solid), the ensem¬ 
bles equally pull each other’s optical phases ^a,b and the light 
phase relaxes to A'ljj = 51(3)°, close to (A'i/jg, = 

(56°, 55°). In contrast, in the more imbalanced (open) data, 
the unobserved (jja is pulled more rapdily toward the phase 
of (ph, and the phase relaxes toward A’ljj — 71(2)°, while 
(A^/jg, Aijjn) = (79°, 73°), i.e., closer to the phase of ensem¬ 
ble b at Tevoi = 0. 

Synchronization necessarily implies moving from a state of 
higher entropy to a state of lower entropy, requiring dissipa¬ 
tion into a bath of states that absorb the entropy. In our atom- 
cavity system, the dominant dissipation mechanism for syn¬ 
chronization is the spontaneously scattered optical pumping 
light involved in re-exciting the atoms from |),) to |t) at rate 
W. Because our atomic ensemble is optically thin in the di¬ 
rection transverse to the laser cavity, the scattering process for 
the ith atom is not collective and causes single-atom collapse, 
erasing the relative quantum phase pi in the single atom su¬ 
perposition state; cos(0i/2) |ti) -f e*'^’ sin(0i/2) —)• |ti). 

It is this relative phase pi that encodes the phase of the single¬ 
atom dipole and thus the phase of the light pi = pi + co'nst 
that is radiated by the single oscillator. It is helpful to vi¬ 
sualize pi as the azimuthal phase of the single-diople Bloch 
sphere and the angle 9i as a polar angle. Most importantly, 
the quantum collapse serves to erase any relative phase error 
Api — pi — (/)avg that had accumulated in time between the 
individual atom’s optical dipole and an appropriately defined 
average of the phases of all of the optical dipoles of partici¬ 
pating atoms ())avg. 

The total cavity field is the sum of the optical fields radiated 
by each atom, with a resulting phase p^wg = pa.vg- This cav¬ 
ity field aligns the optical dipole phase of a newly repumped 
atom to pavg. The combination of realignment to the aver¬ 
age, accrual of phase errors, and erasure of phase errors is the 
physical origin of the quantum synchronization process. 

Another dissipation channel in the atom-cavity system is 
the emission of photons from the cavity mode through the 
mirrors. However, this channel only provides collective infor¬ 
mation to the environment and should not erase single-atom 
phase errors. Detection of a photon exiting the cavity indi¬ 
cates that one atom has made a transition from |t) to j),), but 
it does not indicate which atom made the transition. Still, the 
fast dissipation (k ^ 7j^) ensures that the cavity field adia- 
batically follows the total optical dipole moment of the atoms 
such that the phase of the optical dipoles p^vg can be deter¬ 
mined directly from the phase of the emitted light V’avg- 

We next consider the case in which a continuous source of 
phase error is introduced between the two ensembles. This 
is equivalent to detuning the dressing laser frequencies. As 
S deviates from zero, the total power emitted by the two en¬ 
sembles decreases as shown in Fig. For |i5| > W the to¬ 
tal output power is roughly constant. At the transition point, 
the two ensembles largely behave independently, emitting at 
their respective natural lasing frequencies. The characteristic 


frequency scale is set by W since any relative phase accu¬ 
mulated between the ensembles is reset by repumping. The 
observed maximum synchronized power output is a factor of 
2.2(1) greater than the unsynchronized power output, while 
we predict a factor of 1.8(2). This estimate is based on the 
quenching behavior of the output power with repumping rate 
that accounts for changes in population inversion of each en¬ 
semble II23II27I . The asymmetry of the total power for positive 
and negative 6 is also reflected in the asymmetric behavior in 
the spectra of Fig. and discussed below. 



FIG. 3. (Color online) Total power output vs. detuning for the data 
shown in Fig.Qa). Vertical dashed lines are at the repumping rate 
±lT/27r. Horizontal dashed line is the predicted maximum synchro¬ 
nized output power. 

We can observe the transition from synchronized to unsyn¬ 
chronized behavior in the frequency domain by making het¬ 
erodyne measurements of the light emitted from the cavity. In 
the spectrograms of Fig. each row is a frequency spectrum 
of emitted light from the cavity, with brighter colors indicating 
higher power. Each power spectrum is calculated from 80 /is 
of the time record. The two-dimensional power spectrum is 
created by repeating the measurement at a series of different 
detunings 6, with values shown along the vertical axis. 

For |(5| 3> W, the two ensembles of atoms emit at frequen¬ 
cies very close to the unperturbed Raman transition frequen¬ 
cies. As |5| decreases, the emission frequencies are pulled 
toward each other as the rate of relative phase error introduc¬ 
tion 6 nears the error erasure rate W. We note that we do not 
observe nor expect a region of repulsive synchronization that 
appears when injection locking a single superradiant ensem¬ 
ble to an externally applied drive ll28l . For |(5| < W, the era¬ 
sure of phase errors dominates and the two ensembles radiate 
at a single frequency. 

The observed spectrum qualitatively agrees with the same 
mean-field model used to describe the time dynamics of resyn¬ 
chronization, exhibiting a hyperbolic-like approach (Region 
I) to the synchronized state (Region II). However, there is 
significant asymmetry in the power spectrum. Part of this 
asymmetry arises from a finite detuning of the average Ra¬ 
man transition frequency from resonance with the cavity res¬ 
onance frequency by amount Sc = —2 tt x 4 MHz « k/ 3, an 
operating condition favorable for suppressing relaxation os¬ 
cillations ll26ll^ yet one that introduces an imbalance in the 
coupling to the cavity between ensembles. Other causes of 
asymmetry are imbalances in the optical dipole magnitudes 
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FIG. 4. (Color) Spectrograms of light emitted from two superradiant 
ensembles. Vertical axis is the detuning of dressing lasers 5 and hor¬ 
izontal axis is the Fourier frequency of each power spectrum. The 
power (color scale) is normalized to the maximum power across the 
entire spectrogram, (a) Each power spectrum displayed here repre¬ 
sents the mean of 5 power spectra at each <5. Collective dipoles are 
roughly balanced with Na/Nf, = 0.6 and 7a /ts = 0.8 (b) Asym¬ 
metric operating conditions, Na/Ni, = 1.1 and 'fa/'yb = 1.6. 


(both N and 7) for the data in Fig. Numerical modeling 
indicates that the effects of these small asymmetries are mag¬ 
nified by the interaction between the ensembles. 

We also show in Fig. I^b) that many different behaviors 
can be observed depending on the operating parameters. This 
data shows a significant asymmetry in the emitted power (iii) 
from each ensemble for <5 > 0 and 5 < 0. Many of these 
behaviors are observed in numerical mean-field models of our 
system, but other features, indicated in Fig. [^b) are not: (i) 
the parallel-running frequency component in the lower right 
hand quadrant, (ii) the extra frequency components at ±i5/27r, 
and the asymmetry in the observed linewidth of the two emis¬ 
sion peaks of both Fig. j^a) and (b). The fractional power in 
each sideband (ii) is small, < 8% of the total power in each 
spectrum. 

In prior studies, linewidth broadening was seen to arise 
from an inversion-dependent frequency-pulling mechanism 
that here would cause a common broadening of both 
peaks ||23l|26l. Attempts to identify other classical mecha¬ 


nisms for the asymmetric broadening have been unsuccessful 
and the broadening phenomenon remains an interesting topic 
for future theoretical and experimental study, with the intrigu¬ 
ing possibility that this is a fundamental quantum noise ef¬ 
fect m. 

This work emphasizes the key role that dissipation via re¬ 
pumping plays in erasing the phase errors that are gener¬ 
ated during the synchronization process. This work may 
apply to future technologically relevant implementations of 
superradiant ensembles that would produce optically narrow 
light 1^ . For instance, two superradiant lasers that oper¬ 
ate on atomic transitions with opposite sensitivity to magnetic 
fields could use synchronization to cancel frequency transition 
noise as must be done sequentially in passive optical lattice 
clocks OOl . Also, this work points toward a coupled atom- 
cavity system for exploring quantum noise in phase transi¬ 
tion models m, using atomic synchronization for enhanced 
Ramsey spectroscopy ED, and overcoming the Dick effect by 
transferring coherence between atomic clocks ll^ . 
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